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We numerically investigate the impact on the two-body range of several Newtonian 
and non-Newtonian dynamical effects for some Earth-planet (Mercury, Venus, Mars, 
Jupiter, Saturn) pairs in view of the expected cm-level accuracy in some future planned 
or proposed interplanetary ranging operations. The general relativistic gravitomagnetic 
Lense-Thirring effect should be modeled and solvcd-for in future accurate ranging tests 
of Newtonian and post-Newtonian gravity because it falls within their mcasurability do- 
main. It could a-priori "imprint" the determination of some of the target parameters 
of the tests considered. Moreover, the ring of the minor asteroids, Ceres, Pallas, Vesta 
(and also many other asteroids if Mars is considered) and the Trans-Neptunian Objects 
(TNOs) act as sources of nonncgligible systematic uncertainty on the larger gravitoelec- 
tric post-Newtonian signals from which it is intended to determine the parameters 7 
and 8 of the Parameterized Post-Newtonian (PPN) formalism with very high precision 
(several orders of magnitude better than the current 10 — 4 — 10~ 5 levels). Also other 
putative, nonconvcntional gravitational effects like a violation of the Strong Equivalence 
Principle (SEP), a secular variation of the Newtonian constant of gravitation G, and the 
Pioneer anomaly are considered. The presence of a hypothetical, distant planetary-sized 
body X could be detectable with future high-accuracy planetary ranging. Our analysis 
can, in principle, be extended also to future interplanetary ranging scenarios in which 
one or more spacecraft in heliocentric orbits are involved. The impact of fitting the initial 
conditions, and of the noise in the observations, on the actual detectability of the dy- 
namical signatures investigated, which may be partly absorbed in the estimation process, 
should be quantitatively addressed in further studies. 
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1. Introduction 

Recent years have seen increasing efforts towards the implementation of the Plane- 
tary Laser Ranging (PLR) technique accurate to cm-level 1 ' 2,3 ' 4,5 ' 6,7 ' 8 . It would allow 
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to reach major improvements in three related hclds: solar system dynamics, tests of 
general relativity and alternative theories of gravity, and physical properties of the 
target planet itself. In principle, any solar system body endowed with a solid surface 
and a transparent atmosphere would be a suitable platform for a PLR system, but 
some targets are more accessible than others. Major efforts have been practically 
devoted so far to Mercury 1 and Mars 2 ' 5 , although simulations reaching 93 a.u. or 
more have been undertaken as well 4 ' 7 . In 2005 two interplanetary laser transponder 
experiments were successfully demonstrated by the Goddard Geophysical Astro- 
nomical Observatory (GGAO). The first utilized the nonoptimized Mercury Laser 
Altimeter (MLA) on the Messenger spacecraft 1 ' 3 , obtaining a formal error in the 
laser range solution of 0.2 m, or one part in 10 11 . The second utilized the Mars 
Orbiting Laser Altimeter (MOLA) on the Mars Global Surveyor spacecraft 9,3 . A 
precise measure of the Earth-Mars distance, measured between their centers of mass 
and taken over an extended period (five years or more), would support, among other 
things, a better determination of several parameters of the solar system. Sensitivity 
analyses point towards measurement uncertainties between 2 1 mm and 100 mm. 
The perspectives in measuring the distance between an Earth-based station and an 
active laser transponder on the Martian moon Phobos capable of reaching mm-level 
range resolution have recently been investigated in Ref. 10. The authors of Ref. 11 
envisage the possibility of using also Deimos, in addition to Phobos, as a target 
for an in-situ lander in the framework of the proposed Gravity Experiment with 
TimE Metrology on Martian satEllites (GETEMME) mission a . Its goal is to use 
laser to measure intermartian distances with an accuracy of a few tenth of mm, 
and its nominal duration should be 3 — 5 yr. See also Ref. 12 for an earlier, prelim- 
inary study on the possibility of using both artificial and natural satellites to test 
general relativity in the martian system. Concerning Mercury, a recent analysis on 
the future BepiColombo b mission, aimed to accurately determining, among other 
things, several key parameters of post-Newtonian gravity and the solar quadrupole 
moment from Earth-Mercury distance data collected with a multi-frequency ra- 
dio link 13 ' 14 , points toward a maximum uncertainty of 4.5 — 10 cm in determining 
the Earth-Mercury range over a multi-year time span 13 ' 15,14 (1-8 yr). A proposed 
spacecraft-based mission aimed to accurately measure also the general relativistic 
gravitomagnetic field of the Sun and its adimensional quadrupole mass moment J2 
along with other PPN parameters like 7 and (3 by means of interplanetary ranging is 
the Astrodynamical Space Test of Relativity using Optical Devices (ASTROD) 17 . 
Another space-based missions proposed to accurately test several aspects of the 
gravitational interaction via interplanetary laser ranging are the Laser Astrometric 



a See on the WEB http://meetingorganizer.copernicus.org/EPSC2010/EPSC2010-60.pdf. 
b It is an ESA mission, including two spacecraft, one of which provided by Japan, to be put into 
orbit around Mercury. The launch is scheduled for 2014. The construction of the instruments is 
currently ongoing. 

c Its cheaper version ASTROD I makes use of one spacecraft in a Venus-gravity-assistcd solar orbit, 
ranging optically with ground stations 16 . 
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Test of Relativity (LATOR) 18 , and the interplanetary range to Phobos 10 one goal 
of which is a measurement of 7 with an accuracy of 10 -7 . For a review of the 
motivations for accurately determining the parameters of post-Newtonian gravity, 
in particular (3 and 7, see, e.g., Ref. 19, 20 and references therein. 

In this paper we study the effects that several Newtonian and non-Newtonian 
dynamical features of motion have on the two-body range for the Earth and some 
planets of the solar system for which accurate ranging to spacecraft exists or is 
planned in future. Our goal is to inspect the potential aliasing posed by other com- 
peting dynamical forces acting as source of systematic uncertainty. Indeed, it must 
be recalled that in the range observables actually used in testing post-Newtonian 
gravity there is also a part due to the Earth-planet orbital motions in addition to 
the purely post-Newtonian Shapiro delay connected with the propagation of elec- 
tromagnetic waves. Thus, reaching unprecedented accuracy in only measuring the 
latter effect is useless if the accuracy of the orbital signal is not at a comparable 
level. On the other hand, it should be remarked that we do not aim to quantitatively 
assess the actual measurability of the dynamical effects investigated. It is a different 
and important task which would deserve a dedicated work. Indeed, a fit of the initial 
conditions to the real observations would be needed in order to realistically evaluate 
the level of removal of the effects of interest from the signatures. It is a nontrivial 
task which is beyond the scopes of the present analysis which could be fruitfully 
used to single out the most relevant dynamical signals and focus future efforts on 
them. Anyway, we will use in the following a heuristic rule-of-thumb. The paper 
is organized as follows. In Section 2 we outline the strategy followed and mention 
the Newtonian and non-Newtonian effects investigated. In Section 3, Section 4, Sec- 
tion 5, Section 6, and Section 7 we deal with the ranges of Mercury, Venus, Mars, 
Jupiter and Saturn, respectively. Section 8 is devoted to the conclusions. 

2. Method 

In order to numerically obtain the effect of a given gravitational acceleration, con- 
sidered as a relatively small perturbation P of the Newtonian Sun's monopole, on 
the range \p\ between the Earth-Moon Barycenter (EMB) and a planet we used 
MATHEMATICA to simultaneously integrate with the Runge-Kutta method the 
equations of motion in Cartesian coordinates of EMB and the planet considered 
with and without the perturbation P investigated by using the same set of initial 
conditions. We adopted the ICRF/J2000.0 reference frame, with the ecliptic and 
mean equinox of the reference epoch, i.e. J2000.0, centered at the Solar System 
Barycenter (SSB); the initial conditions at the epoch J2000.0 were retrieved with 
the HORIZONS WEB interface by JPL, NASA. The temporal interval of the numer- 
ical integration for Mercury and Venus has been taken equal to At = 2 yr in view of 
the fact that the typical operational time spans envisaged for future PLR technique 
are similar. For Mars and Saturn, for which ranging to spacecraft is currently on- 
going, we adopted At = 5 yr. Also for Jupiter we used At — 5 yr. The basic model 
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adopted consists of the barycentric equations of motion of the Sun, the eight plan- 
ets, the Moon, Ceres, Pallas, Vesta, Pluto and Eris, to be simultaneously integrated; 
the forces acting on them include the mutual Newtonian TV— body interactions, the 
perturbation due to the solar quadrupolar mass moment J2, the effect of two rings 
modeling the actions of the minor asteroids and of the Trans-Neptunian Objects 
(TNOs), and the general relativistic gravitoelectric Schwarzschild and gravitomag- 
netic Lense-Thirring fields of the Sun. As additional perturbations, we modeled the 
action of a distant, planet-like body kept fixed in a given spatial position (planet 
X), secular rate G/G of the Newtonian gravitational constant G, a violation of 
the Strong Equivalence Principle (SEP) by means of the Nordtvedt parameter 77, a 
constant and uniform acceleration radially directed towards the Sun, and acting on 
Uranus, Neptune, Pluto (and Eris), to account for the Pioneer anomaly. 

In order to preliminarily give an idea of the potential measurability of the effects 
considered, the computed differences A|p| = \pp\ — |/?r|, where R refers to a reference 
orbit which does not contain the perturbation P of interest, were subsequently 
compared to the available time series of the range residuals for the inner planets and 
Saturn which set the present-day accuracy level in ranging to planets 21,22 . When the 
possibility that a given, unmodeled dynamical effect may show or not its signature 
in the range residuals it must be considered that the magnitude of such an effect 
should roughly be one order of magnitude larger than the range residuals accuracy. 
This to avoid the risk that it may be absorbed and partially or totally removed from 
the signature in the process of estimation of the initial conditions and of the other 
numerous solve-for parameters in the real data reduction. As outlined before, it is 
just a sort of rule of thumb; fully quantitative, realistic analyses, outside the scopes 
of the present paper, should require the actual fitting of the initial conditions to the 
observations. 

Depending on the dynamical effect one is interested in, some of the perturbations 
examined here are to be considered as sources of noise inducing systematic bias on 
the target signal. For example, if the goal of the analysis is, say, the Lense-Thirring 
effect, then the range perturbation due to the TNOs is clearly a source of potential 
systematic error which has to be evaluated. Thus, our plots are useful to assess 
the level of aliasing of several potential sources of aliasing for some non-Newtonian 
effects and the correlations that may occur in estimating them. Dynamical effects 
which are viewed as noise in a given context can also be regarded as main targets in 
another one; see, e.g., the proposed determination of asteroid masses through the 
ASTROD mission 23 . 

3. Earth-Mercury range 

At present, the 1-way range residuals of Mercury from radar-ranging span 30 yr 
(1967-1997) and are at a few km-level (Figure B-2 of Rcf. 21); the same holds for 
the 1-way Mercury radar closure residuals covering 8 yr (1989-1997, Figure B-3 a) of 
Ref. 21). There are also a pair of Mariner 10 range residuals in the 70s at Mercury at 
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0.2 km level (Figure B-3 b) of Ref. 21). Ranging to BcpiColombo should be accurate 
to 15,14 4.5 — 10 cm over a few years. 



3.1. The Schwarzschild field of the Sun 

In Figure 1 we plot the effect of the gravitoelectric Schwarzschild field of the Sun 
on the Earth-Mercury range. We modeled its acceleration as 24 

r + 2(1 + 7) (r-i>)t>|, (1) 

where we inserted the PPN parameters (3 and 7: they are equal to 1 in general 
relativity and we used such values. Figure 1 can be compared with Figure 1 of 
Ref. 14, obtained for unspecified initial conditions' 1 : they arc quite similar. The 
maximum variation of the signal is of the order of 4 x 10 5 m, corresponding to 
a measurement accuracy of about 2.5 x 10 -7 . The expected realistic accuracy in 
determining f3 and 7 is 2 x 10 -6 in BepiColombo 13 . 



ass = 



GM 



GM 

2(p + 7) jv • v 



3.2. The oblateness of the Sun 

Figure 2 shows the nominal effect of the Sun's quadrupolar mass moment on the 
Mercury range for J2 = 2 x 10~ 7 . Its action has been modeled as 25 

r + 2(f-fc)fc|, (2) 

where R is the Sun's mean equatorial radius and k is the unit vector of the z axis 
directed along the body's rotation axis. Since eq. (2) holds in a frame with its {xy} 
plane coinciding with the body's equator, we rotated the the mean ecliptic at the 
epoch to the Sun's equator which is inclined to it by the Carrington angle 26 i = 7.15 
deg. The signal of Figure 2 has a maximum span of 300 m, corresponding to an ac- 
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6J 2 n UM 

2^ 



i l-5(r-fc) 2 



d It also includes the Shapiro delay contribution. 
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Fig. 1. Difference A\p\ = |pp| — |pr| in the numerically integrated EMB-Mercury ranges with 
and without the perturbation due to the Sun's Schwarzschild field over At = 2 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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curacy measurement of 3 x 10~ 4 . A measure of the solar J2 accurate to 10~ 2 is one 
of the goals of BepiColombo 13 ; knowing precisely J2 would yield important insights 
on the internal rotation of the Sun. At present, it is known with an uncertainty of 
about 29 10%. The solar quadrupole mass moment may play the role of source of 
systematic bias with respect to, e.g., some non- Newtonian dynamical effects. Con- 
cerning the gravitoelectric signal previously analyzed, the mismodeled J2 signature 
would impact it a 7.5 x 10~ 5 level. It is important to note that the patterns of the 
two signals are rather different. Conversely, as we will see, the determination of J 2 
at the desired level of accuracy may be affected by other unmodeled/mismdeled 
dynamical effects acting as systematic sources of aliasing on it. 



3.3. The Lens e-Thir ring effect of the Sun 

Figure 3 depicts the range perturbation due to the Sun's Lense-Thirring effect, nei- 
ther considered so far in the dynamical force models of the planetary ephemerides 
nor in the BepiColombo analyses. It is a general relativistic feature of motion in- 
duced by the rotation of the Sun which acts upon a test particle moving with 
velocity v with a noncentral acceleration 24 



(1+ 7 )G 



3 

— (r x v) (r • S) + (v X S) 



(3) 



where S is the Sun's proper angular momentum. According to hclioscismology ,28 , 
its magnitude is S — (190.0 ± 1.5) x 10 39 kg m 2 s™ 1 . Note that such a value does 
not come from planetary orbital dynamics, so that there is no risk of a-priori "im- 
printing" of general relativity itself on range tests of the solar Lense-Thirring effect 
which could, thus, be regarded as genuine and unbiased. Also in this case we ro- 
tated the reference frame to the mean ecliptic at the epoch to the Sun's equator 
by the Carrington angle because eq. (3) holds in a frame with its z axis aligned 
with S. The peak-to-peak amplitude of the Lense-Thirring signal is up to 17.5 m 

EMB-Mercury rLui^e: SSB iiuiuericalcalculation 




Fig. 2. Difference A|p| = |pp — \pr\ in the numerically integrated EMB-Mercury ranges with and 
without the nominal perturbation due to the Sun's quadrupole mass moment J2 = 2.0 X 10 — 7 over 
At = 2 yr. The same initial conditions (J2000.0) have been used for both the integrations. The 
state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. The 
integrations have been performed in the ICRF/J2000.0 reference frame, with the mean equinox of 
the reference epoch and the reference {xy} plane rotated from the mean ecliptic of the epoch to 
the Sun's equator, centered at the Solar System Barycenter (SSB). 
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over 2 yr, which, if on the one hand is unmcasurable from currently available radar- 
ranging to Mercury, on the other hand corresponds to a potential relative accuracy 
in measuring it with BepiColombo of 2 — 5 x 10 -3 ; this clearly shows that the solar 
gravitomagnctic field should be taken into account in future analyses and data pro- 
cessing. Otherwise, it would alias the recovery of other effects. For example, it may 
affect the determination of J2 at 12% level. On the other hand, in order to allow for 
a determination of the Lense-Thirring effect, the Sun's quadrupole mass moment 
should be known with an accuracy better than the present-day one by at least one 
order of magnitude; this is just one of the goals of BepiColombo. Moreover, since 
the Lense-Thirring effect depends on c 7, neglecting it may alias the determination 
of 7 through the larger gravitoelectric signal at 4 x 10 -5 level. 



3.4. The ring of the minor asteroids and Ceres, Pallas and Vesta 

In Figure 4 we depict one potential source of systematic bias, i.e. the action of the 
ring of minor asteroids 29 . We modeled it following Ref. 30. For those planets for 
which r > R lmg , by posing a = R T i ns /r, we obtained 

Gr77l r j n g ( n 3 9 45 



from 



2r 3 V 2 32 
Gm r i 



H-a 2 + -a Mr, (4) 



-inner ring 2r 3 

For r < i? r ing, by posing a = r/R r - mg , we obtained 

Gm r ing / 9 3 75 5 



ab^(a)-bf(a) r. (5) 



"ring 



Oouter ring 2rR2 I a + -a + — a ) r, (6) 



c The multiplicative factor 2 in front of eq. (3) comes from 24 1 + 7. 
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Fig. 3. Difference A|p| = |pp| — |pr | in the numerically integrated EMB-Mercury ranges with and 
without the perturbation due to the Sun's Lense-Thirring field over At = 2 yr. The same initial 
conditions (J2000.0) have been used for both the integrations. The state vectors at the reference 
epoch have been retrieved from the NASA JPL Horizons system. The integrations have been 
performed in the ICRF/J2000.0 reference frame, with the mean equinox of the reference epoch 
and the reference {xy} plane rotated from the mean ecliptic of the epoch to the Sun's equator, 
centered at the Solar System Barycenter (SSB). 
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ring 



r. 



(7) 



Recall that the Laplace coefficients are defined as 

1 [ 2lT cos jtpdtp 



(1 — 2a cos?/; + a 2 ) s 1 
where s is a half-integer; a useful approximate expression in terms of a series is 31 



(8) 



s(s + l)...(s + j- 1) 
1-3---J 



-or 



;i+i) 



s(s + l)(s 



a 



l-2(j + l)(j + 2) 



(9) 



By assuming for the ring of the minor asteroids a nominal mass of 29 m r j ng = 
1 x 10~ 10 M Q and a radius 29 i? r i ng = 3.14 au, it would impact the Mercury range at 
4 m level (peak-to-peak amplitude), which is, in fact, measurable. Its nominal bias on 
the Schwarzschild, J2 and Lense-Thirring signals would be 1 x 10~ 5 , 1.3 x 10 -2 , 2.3 x 
10 _1 , respectively. Anyway, the present-day level of uncertainty in the mass of the 



ring is <5m r 



0.3 x 10 M©. Thus, the impact of such a mismodeling would 



be, 3 x 10~ 6 ,4 x 10 _3 ,7 x 10 -2 , respectively; it cannot be considered negligible. 

The effect of Ceres, Pallas and Vesta on the determination of some Newtonian 
and non-Newtonian parameters with BepiColombo has been preliminarily inves- 
tigated in Rcf. 15. Here in Figure 5 we show the nominal perturbation on the 
Earth-Mercury range due to the combined actions of Ceres, Pallas and Vesta; the 
values for their masses have been retrieved from Ref. 32. Its peak-to-peak amplitude 
amounts to 80 m; thus, their signature would be measurable at a 0.6 — 1 x 10~ 3 
level. Anyway, the mismodeled solar quadrupole mass moment would bias their 
signal at 4 x 10 _1 level. The Lense-Thirring effect, if unmodeled, would have an 
impact at 2.2 x 10 _1 level. The present-day relative uncertainties in their masses 
are 6 x 10 _3 ,3 x 10~ 2 ,2 x 10~ 2 respectively 32 . This implies a mismodeled signal 



EMB-Mercury range SSB iiumericalcalculaiion 




Fig. 4. Difference A|p| = |pp| — \pn\ in the numerically integrated EMB-Mercury ranges with and 
without the nominal perturbation due to the ring of minor asteroids with 29 m r i ng = 1 x 10 — '"Mq 
and -R r i n g = 3.14 a.u. over At = 2 yr. The same initial conditions (J2000.0) have been used for 
both the integrations. The state vectors at the reference epoch have been retrieved from the NASA 
JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 reference frame, 
with the ecliptic and mean equinox of the reference epoch, centered at the Solar System Baryccnter 
(SSB). 
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with a peak-to-peak amplitude of 50 cm. It would impact the Schwarzschild, Ji and 
Lense-Thirring range perturbations at 1 x 10~ 6 , 2 x 10~ 3 , 3 x 10 -2 level, respectively. 

3.5. The Trans- Neptunian Objects 

The situation is different for another potential source of systematic uncertainty, i.e. 
the Trans- Neptunian Objects (TNOs). Figure 6, obtained by modeling them as a 
ring with 22 m r i n g = 5.26 x 10 -8 M Q and 22 -R r i ng = 43 au, shows that their maximum 
effect would amount to 80 cm. We used the same formulas as for the asteroid ring. 
Such an effect, not taken into account so far, would be better measurable than that 
by the minor asteroids. This implies a bias of 2 x 10~ 6 on the Schwarzschild signal, 
3 x 10 -3 for J2 and 4.5 x 10 -2 for the Lense-Thirring effect. A major concern is that 
the mass of the TNOs is far from being accurately known, so that an uncertainty 
as large as 100% should be applied. 



EMB-Mercury range: SSB numerical calculation 
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Fig. 5. Difference A\p\ = |pp| — |ppj in the numerically integrated EMB-Mercury ranges with 
and without the nominal perturbation due to 32 Ceres, Pallas, Vesta over At = 2 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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Fig. 6. Difference A|p| = |pp| — |pr in the numerically integrated EMB-Mercury ranges with 
and without the nominal perturbation due to the ring of Trans-Neptunian Objects with 22 m r ; ng = 
5.26 X 10 — 8 Mq and firing = 43 a.u. over At = 2 yr. The same initial conditions (J2000.0) have 
been used for both the integrations. The state vectors at the reference epoch have been retrieved 
from the NASA JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 
reference frame, with the ecliptic and mean equinox of the reference epoch, centered at the Solar 
System Barycenter (SSB). 
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3.6. Violation of the Strong Equivalence Principle 

Let us, now, focus our attention to other nonstandard effects like a SEP violation, 
a variation of the Newtonian gravitational constant, the Pioneer anomaly and a 
putative planet X. 

Concerning the SEP violation, we modeled it in the acceleration a, of a body i 

as 



(SEP) 



r (g) - 




(i) 
_ m\ _ 


SEP \ 



Grrij 



1,2,. ..N. 



(10) 



In it and 
body i. Their ratio is 



are the gravitational and inertial masses, respectively, of the 



1 + TjQj 



(11) 



J SEP 



in which r\ is the Nordtvedt dimensionless constant 33,34 accounting for SEP viola- 
tion, and 



a 



rriiC 2 



(12) 



where Ei is the (negative) gravitational self-energy of the i— th body, and m^c 2 is 
its total mass-energy. For a spherical body of radius Ri 19 , 



n 3 Gm, 

n< = -■=- 



5 Ric 2 ' 



(13) 



Figure 7 shows its effect on the Mercury range for rj = 10~ 5 ; at present, the most 
accurate constraints on r\ come from Lunar Laser Ranging (LLR) amounting to 35 
rj = (4.0 ± 4.3) x 10 -4 . The peak-to-peak amplitude of the signature of Figure 7 is 
6 mm, which is practically impossible to detect. Moreover, the SEP signal would 
be totally swamped by other dynamical effects like the ones by the minor asteroid 
ring, even if modeled at the present-day level of accuracy, and the TNOs. Also the 
Lensc-Thirring effect, modeled or not, would be of concern. 



3.7. Secular variation of the Newtonian constant of gravitation 

The case of a possible variation of the Newtonian gravitational constant is inter- 
esting because recently E.V. Pitjcva in Rcf. 22 preliminarily reported a secular 
variation for it G/G = (—5.9 ± 4.4) x 10 -14 yr _1 , statistically significant at 3 — a 
level. Other researchers get results statistically compatible with 0. Folkner in Ref. 36 
gets an upper bound of 2 x 10~ 13 yr _1 from planetary ephemerides as well. Williams 



In terms of the PPN parameters fi and 7 it is r\ = A/3 — 7 — 3, so that rj = in general relativity. 
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et al. in Ref. 35 obtain G/G = (4±9) x 10 



-13 



yr 



from LLR. We modeled a secular 



variation of G in the equations of motion according to 

G [l + g(t-£ ) 



a 



(G) 



E 



i J ■ 



i = 1, 2, ...N. 



(14) 



Figure 8 depicts its impact on the Mercury range. Its maximum effect would be 
about 60 cm, which, in principle, should be measurable with BcpiColombo at a 
0.7 — 2 x 10 _1 level of relative accuracy. If a secular decrease of G will be confirmed as 
a genuine physical effect by further analyses of planetary data by independent teams 
of astronomers, it should be modeled in the BcpiColombo data analysis because, 
otherwise, it would affect the Schwarzschild, J2 and Lense-Thirring signatures at 
1 x 10~ 6 ,2 x 10~ 3 ,3 x 10~ 2 level, respectively. Anyway, the action of the Lense- 
Thirring effect, of Ceres, Pallas, Vesta and of TNOs, modeled or not, would likely 
bias the recovery of the putative G signal in a severe way; it must be recalled that 
the mismodeled signature due to the three large asteroids is 50 cm. 
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Fig. 7. Difference A|p| = |pp| — IprI in the numerically integrated EMB-Mercury ranges with and 
without the perturbation due to a violation of the Strong Equivalence Principle for r\ = 10" 5 over 
At = 2 yr. The same initial conditions (J2000.0) have been used for both the integrations. The 
state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. 
The integrations have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and 
mean equinox of the reference epoch, centered at the Solar System Barycenter (SSB). 
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Fig. 8. Difference A|p| = \pp\ — |pr,| in the numerically integrated EMB-Mercury ranges with and 
without the nominal perturbation due to a secular variation of G as large as 22 G/G = —5.9 X 10" 14 
yr -1 over At = 2 yr. The same initial conditions (J2000.0) have been used for both the integrations. 
The state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. 
The integrations have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and 
mean equinox of the reference epoch, centered at the Solar System Barycenter (SSB). 
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3.8. The Pioneer Anomaly 

The Pioneer anomaly 37 , which is a constant anomalous extra-acceleration approxi- 
mately directed towards the Sun of magnitude Api = 8.74 x 10~ 10 m s -2 detected 
in the telemetry of the Pioneer 10/11 spacecraft after they passed 20 au, may, in 
principle, impact the Earth- Mercury range as well in an indirect way through the al- 
tered action on them of the bodies directly affected by such a putative exotic force, 
i.e. Uranus, Neptune, Pluto and Eris. This would be another way of testing the 
hypothesis of a gravitational nature of the Pioneer anomaly in addition to directly 
looking at the outer planets 38 which has given negative results 39 ' 29 . In Figure 9 we 
plot its signature. It would amount to 4 mm, which is too small to be realistically 
detected. Moreover, also the aliasing bias of the other effects previously considered 
would be crucial. 

3.9. Planet X 

Finally, let us consider the potential ability of BepiColombo of detecting the sig- 
nature of a putative remote planet X. In Figure 10 and Figure 11 we depict the 
X's range signal for the minimum and maximum value of its tidal parameter 
/Cx = GMx/r| according to the anomalous perihelion precession of Saturn an- 
alyzed in Ref. 40. Let us recall that recently Pitjeva in Ref. 41 and Fienga et al. in 
Ref. 29 independently determined statistically significant extra-precessions of the 
perihelion of Saturn from preliminary analysis of some years of radio-tracking data 
of the Cassini spacecraft. Their values are 8 Arijs a t = —6 ± 2 mas cty -1 , cited in 
Ref. 29, and 29 Arosat = —10 ±8 mas cty" 1 . Further data analyses of longer Cassini 
data records are required to confirm or disproof the existence of such an anomaly as 



g E.V. Pitjeva, private communication, December 2008. 
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Fig. 9. Difference A|p| = |pp — |pp | in the numerically integrated EMB-Mercury ranges with and 
without the nominal perturbation due to a Pioneer-like constant and unform radial acceleration 
of 37 Api = 8.74 X 10~ 10 m s -2 acting upon Uranus, Neptune, Pluto, Eris over At = 2 yr. The 
same initial conditions (J2000.0) have been used for both the integrations. The state vectors at 
the reference epoch have been retrieved from the NASA JPL Horizons system. The integrations 
have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of 
the reference epoch, centered at the Solar System Barycenter (SSB). 
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a genuine physical effect' 1 . The maximum effect of X on the Mercury range would 
be as large as 1.5 — 3 m. falling within the measurability domain of BcpiColombo. 
It must be noted that a large part of such a signal may be largely confused with the 
action of the TNOs. The Lense-Thirring effect, if not modeled, would be another 
source of serious systematic error as well. 

4. Earth- Venus range 

Although, at present, no tests of interplanetary ranging to Venus have been practi- 
cally performed, contrary to Mercury and Mars, we prefer to treat also its case not 
only for completeness but also because simulations of interplanetary transponder 



h Pitjeva in Ref. 22 reports a new value which, instead, is statistically compatible with 0, i.e. 
Arosat = -10 ± 15 mas cty - 1 
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Fig. 10. Difference A|p| = |pp| — |ppj in the numerically integrated EMB-Mercury ranges with 
and without the perturbation due to hypothetical remote planet X lying almost in the ecliptic with 
minimum tidal parameter 40 /Cx = 1-6 X 10 -26 s — 2 over At = 2 yr. The same initial conditions 
(J2000.0) have been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference epoch, 
centered at the Solar System Barycenter (SSB). 
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Fig. 11. Difference A|p| = |pp| — |ppj in the numerically integrated EMB-Mercury ranges with 
and without the perturbation due to hypothetical remote planet X lying almost in the ecliptic with 
maximum tidal parameter 40 /C x = 2.7 X 10 -26 s~ 2 over At = 2 yr. The same initial conditions 
(J2000.0) have been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference epoch, 
centered at the Solar System Barycenter (SSB). 
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and laser communications experiments via dual station ranging to SLR satellites 
covering also Venus have been implemented 4,7 . Currently available radar-ranging 
normal points to Venus cover about 33 yr, from 1962 to 1995. The range residuals 
are depicted in Figure B-6 of Rcf. 21; after having been as large as 15 km in the first 
10 yr, they drop below 5 km in the remaining. Figure B-4 shows the range residuals 
to Venus Express at Venus from 2006 to 2008; the are below the 10 m level. 



4.1. The Schwarzschild field of the Sun 

According to Figure 12, the peak-to-peak amplitude of the general relativistic 
Schwarzschild effect is 1.2 x 10 5 m over At — 2 m. A hypothetical, future laser rang- 
ing to, say, a suitably equipped target orbiting Venus accurate to 10 cm would allow 
to measure such a general relativistic signal with a relative accuracy of 8 x 10~ 7 . 
In the following we will discuss the corrupting impact of some potential sources of 
systematic errors. 



4.2. The oblateness of the Sun 

The range perturbation due to the Sun's oblateness is depicted in Figure 13 for the 
nominal value J 2 = 2 x 10~ 7 . Also in this case a barycentric frame rotated to the 
Sun's equator has been adopted. The nominal maximum shift is 40 m, so that a 
measure accurate to2.5xl0 -3 would be possible with a future 10 cm-level ranging 
technique. Viewed as a source of systematic uncertainty, the solar quadrupole mass 
moment would affect the Schwarzschild signal at 3 x 10 -6 level by assuming the 
present-day uncertainty in it, i.e. 10%. The temporal patterns of the two signals are 
quite different. Note that the dynamical action of J2 was modeled in producing the 
Venus Express residuals; thus, a mismodeled signal as large as just 4 m should have 
been left, in agreement with the range residuals. 




Fig. 12. Difference A|p| = |pp| — |pr| in the numerically integrated EMB- Venus ranges with 
and without the perturbation due to the Sun's Schwarzschild field over At = 2 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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4.3. The Lense-Thirring effect of the Sun 

Figure 14 shows the Lense-Thirring perturbation of the Venus range, integrated in 
a frame aligned with the Sun's equator. The peak-to-peak amplitude is 2 m. which 
would be measurable with a future accurate cm-level ranging device with a relative 
accuracy of 2 — 5 x 10~ 2 . The Lense-Thirring signature is still too small to be 
detected nowadays with the current spacecraft ranging. If not modeled, the Sun's 
gravitomagnctic field would impact a determination of J 2 at a 5% level, while the 
Schwarzschild signal would be biased by the Lense-Thirring one at a 3 x 10~ 5 level. 
It must be noted that the two relativistic signals exhibit very similar patterns. The 
present-day 10% uncertainty in the Sun's oblateness would yield a mismodeled signal 
two times larger than the gravitomagnetic one. Anyway, their temporal signatures 
are different, so that it would be possible, in principle, to separate them. 
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Fig. 13. Difference A|p] = |pp| — ppj in the numerically integrated EMB- Venus ranges with and 
without the nominal perturbation due to the Sun's quadrupole mass moment J2 = 2.0 X 10 — 7 over 
At = 2 yr. The same initial conditions (J2000.0) have been used for both the integrations. The 
state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. The 
integrations have been performed in the ICRF/J2000.0 reference frame, with the mean equinox of 
the reference epoch and the reference {xy} plane rotated from the mean ecliptic of the epoch to 
the Sun's equator, centered at the Solar System Barycenter (SSB). 
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Fig. 14. Difference A|p| = |pp| — |ppj in the numerically integrated EMB- Venus ranges with and 
without the perturbation due to the Sun's Lense-Thirring field over At = 2 yr. The same initial 
conditions (J2000.0) have been used for both the integrations. The state vectors at the reference 
epoch have been retrieved from the NASA JPL Horizons system. The integrations have been 
performed in the ICRF/J2000.0 reference frame, with the mean equinox of the reference epoch 
and the reference {xy} plane rotated from the mean ecliptic of the epoch to the Sun's equator, 
centered at the Solar System Barycenter (SSB). 
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4.4. The ring of the minor asteroids and Ceres, Pallas and Vesta 

The impact of the ring of the minor asteroids on the Venus range is depicted in 
Figure 15. With its nominal maximum span of 3 m (peak-to-peak amplitude), also in 
this case such a perturbation would be detectable with a cm-level ranging, and may 
pose some problems to the other signals of interest previously examined. Indeed, its 
mismodeled signature would impact the Lense-Thirring one at 4.5 x 10 _1 level, while 
the bias on J 2 and the Schwarzschild effect is 2.3 x 10~ 2 and 7.5 x 10 -6 , respectively. 
Anyway, the time signatures arc different. Figure 16 shows the nominal perturbation 
on the Venus range by Ceres, Pallas, Vesta. The peak-to-peak amplitude is 175 m, 
measurable at a 6 x 10 -4 level with a ranging device accurate to 10-cm. The aliasing 
effect of the current mismodeling in J2 and in the unmodeled Lense-Thirring effect 
is of the order of 2 x 1CT 2 , 1 x 1CT 2 , respectively. Conversely, it turns out that the 
peak-to-peak amplitude of the mismodeled signature of Ceres, Pallas, Vesta is 1 
m. It would largely alias the Lense-Thirring effect, while the systematic relative 
uncertainty induced on the gravitoelectric and J2 range perturbations would be 
8 x 1CT 6 ,2 x 1CT 2 , respectively. 

4.5. The Trans-Neptunian Objects 

Figure 17 shows the effect of the TNOs on the Venus range. Its peak-to-peak am- 
plitude is about 50 cm: it may be detectable. Its bias on the Schwarzschild, J2 and 
Lense-Thirring signals is 4 x 10~ 6 , 1.2 x 10~ 2 ,2.5 x lO^ 1 , respectively. Concerning 
the aliasing effect on the gravitomagnetic effect, it must be noted that the temporal 
evolution of the two signals is different. This would help in decorrclating them. 

4.6. Violation of the Strong Equivalence Principle 

Figure 18 illustrates the nominal perturbation of the venusian range due to a SEP 
violation with rj = 10 -5 . Such a signal is completely negligible because its peak-to- 
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Fig. 15. Difference A|p| = |pp| — in the numericaiiy integrated EMB- Venus ranges with and 
without the nominal perturbation due to the ring of minor asteroids with 29 m r i ng = 1 X 10 — '"Mq 
and i?, r i n g = 3.14 a.u. over At = 2 yr. The same initial conditions (J2000.0) have been used for 
both the integrations. The state vectors at the reference epoch have been retrieved from the NASA 
JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 reference frame, 
with the ecliptic and mean equinox of the reference epoch, centered at the Solar System Baryccnter 
(SSB). 
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peak amplitude is of just 8 mm. Apart from the fact that it would be undetectable, 
it would be overwhelmed by all the other signatures, modeled or not, previously 
considered. 

4.7. Secular variation of the Newtonian constant of gravitation 

The effect of a secular variation of G as large as G/G = —5.9 x 10~ 14 yr -1 is 
depicted in Figure 19. Its peak-to-peak amplitude is about 7 cm, which is hardly 
detectable even with a cm-lcvcl ranging system. Moreover, such a signature would 
be easily biased by the other dynamical effects considered; for example, recall that 
the mismodclcd effect of Ceres, Pallas, Vesta is as large as 1 m. 

4.8. The Pioneer Anomaly 

Figure 20 illustrates the Venus range perturbation induced by the indirect effect of 
the Pioneer anomaly assumed acting only on the outer planets of the solar system. 
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Fig. 16. Difference A|p] = |pp| — ppj in the numericaiiy integrated EMB- Venus ranges with 
and without the nominal perturbation due to 32 Ceres, Pallas, Vesta over At = 2 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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Fig. 17. Difference A|p| = |pp| — |ppj in the numerically integrated EMB- Venus ranges with and 
without the nominal perturbation due to the ring of Trans-Neptunian Objects with 22 m r ; ng = 
5.26 X 10 — 8 Mq and i? r i ng = 43 a.u. over At = 2 yr. The same initial conditions (J2000.0) have 
been used for both the integrations. The state vectors at the reference epoch have been retrieved 
from the NASA JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 
reference frame, with the ecliptic and mean equinox of the reference epoch, centered at the Solar 
System Barycenter (SSB). 
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Their motions would be altered with respect to the standard case, thus reflecting 
also on the Earth- Venus distance. Its peak-to-peak amplitude amounts to 5 mm, 
about equal to the corresponding effect for Mercury (4 mm). It is negligible because 
it would be undetectable, given the expected cm-level accuracy of future inter- 
planetary ranging devices. Moreover, the much larger aliasing effects of the other 
competing dynamical signatures would completely overwhelm it. 

4.9. Planet X 

More interesting is the situation for a putative planet X. Indeed, the peak-to-peak 
amplitude of its signal, illustrated in Figure 21-Figure 22, is 3 — 5 m. It would be 
detectable in future if a cm-level accuracy in ranging will be achieved. For the mo- 
ment, the X's signature is compatible with the currently available range-residuals. 
Note that its pattern would be similar to that due to the TNOs ring, but its mag- 
nitude would be up to 10 times larger. On the contrary, the Lense-Thirring effect 
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Fig. 18. Difference A|p] = |pp| — ppj in the numericaUy integrated EMB- Venus ranges with 
and without the nominal perturbation due to a violation of the Strong Equivalence Principle for 
r\ = 10~ 5 over At = 2 yr. The same initial conditions (J2000.0) have been used for both the 
integrations. The state vectors at the reference epoch have been retrieved from the NASA JPL 
Horizons system. The integrations have been performed in the ICRF/J2000.0 reference frame, with 
the ecliptic and mean equinox of the reference epoch, centered at the Solar System Baryccnter 
(SSB). 
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Fig. 19. Difference A|p] = |pp| — ppj in the numerically integrated EMB- Venus ranges with and 
without the nominal perturbation due to a secular variation of G as large as 22 G/G = —5.9 X 10~ 14 
yr -1 over At = 2 yr. The same initial conditions (J2000.0) have been used for both the integrations. 
The state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. 
The integrations have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and 
mean equinox of the reference epoch, centered at the Solar System Baryccnter (SSB). 
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has a comparable, m-level size, but a different time signature. 



5. Earth-Mars range 

For Mars we have at our disposal long time series of range residuals accurate to 
about 1 — 10 m-level thanks to several spacecraft (Viking, Mars Pathfinder, Mars 
Global Surveyor, Mars Odyssey, Mars Reconnaissance Orbiter, Mars Express) which 
have orbited, or are still orbiting, the red planet. Figure B-10 of Rcf. 21 depicts the 
1-way range residuals of the Viking Lander at Mars spanning from 1976 to 1982; 
they are at approximately 20 m level. Figure B-ll of Ref. 21 shows the 1-way range 
residuals of several post- Viking spacecraft; they generally cover a few years and are 
accurate to 5 — 10 m. In the following, we will adopt an integration time span At = 5 
yr. 
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Fig. 20. Difference A|p] = pp| — ppj in the numerically integrated EMB- Venus ranges with and 
without the nominal perturbation due to a Pioneer-like constant and unform radial acceleration 
of 37 A Pio = 8.74 X 10 -10 m s~ 2 acting upon Uranus, Neptune, Pluto, Eris over At = 2 yr. The 
same initial conditions (J2000.0) have been used for both the integrations. The state vectors at 
the reference epoch have been retrieved from the NASA JPL Horizons system. The integrations 
have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of 
the reference epoch, centered at the Solar System Barycenter (SSB). 
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Fig. 21. Difference A|p| = |pp| — |ppj in the numerically integrated EMB- Venus ranges with and 
without the perturbation due to hypothetical remote planet X lying almost in the ecliptic with 
minimum tidal parameter 40 /Cx = 1-6 X 10 -26 s — 2 over At = 2 yr. The same initial conditions 
(J2000.0) have been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference epoch, 
centered at the Solar System Barycenter (SSB). 
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5.1. The Schwarzs child field of the Sun 

The general relativistic Schwarzschild perturbation of the Mars range is shown in 
Figure 23. Its peak-to-peak amplitude is 2.5 x 10 5 m. Thus, a 5 — 10 cm-level 
ranging device operating continuously over 5 years would allow a relative accuracy 
in measuring it of 2 — 4 x 10~ 7 . Also in this case, several potential sources of 
systematic errors are to be carefully considered. 

5.2. The oblateness of the Sun 

Figure 24 illustrates the nominal signal due to the Sun's quadrupole mass moment 
for J2 = 2 x 10 -7 computed in a frame aligned with the Sun's equator. Its peak-to- 
peak amplitude amounts to about 70 m; thus, its effect would be well measurable 
at a 7 x 10 -4 — 1 x 10~ 3 level by means of a new ranging facility with an accuracy of 
the order of cm. Concerning its actual presence in the present-day range residuals, 
it must be noted that the dynamical action of the solar J2 has always been mod- 
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Fig. 22. Difference A|p| = (pp| — |ppj in the numerically integrated EMB- Venus ranges with and 
without the perturbation due to hypothetical remote planet X lying almost in the ecliptic with 
maximum tidal parameter 40 /Cx = 2.7 x 10~ 26 s~ 2 over At = 2 yr. The same initial conditions 
(J2000.0) have been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference epoch, 
centered at the Solar System Barycenter (SSB). 
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Fig. 23. Difference A|p| = |pp| — |ppj in the numerically integrated EMB-Mars ranges with 
and without the perturbation due to the Sun's Schwarzschild field over At = 5 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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eled in producing them. Since J2 is nowadays accurate to lO -1 , the corresponding 
mismodeled signature would be as large as about 7 m. i.e. well compatible with 
the range residuals available. Its impact as a source of systematic uncertainty on 
the Schwarzschild signal amounts to 3 x 10 -5 ; note, however, the different time 
signatures of Figure 23 and Figure 24. 

5.3. The Lense-Thirring effect of the Sun 

The Lense-Thirring range perturbation, computed in a frame aligned with the Sun's 
equator, is shown in Figure 25. Its peak-to-peak amplitude is about 4 m, not too 
far from the present-day range accuracy; thus, its existence as predicted by general 
relativity is not in contrast with the range residuals currently available. It could 
be measured with a future cm-level ranging system at a 1 — 2.5% level. If not 
properly modeled, the gravitomagnetic signature would impact the Schwarzschild 
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Fig. 24. Difference A|p| = \pp\ — |pr| in the numerically integrated EMB-Mars ranges with and 
without the nominal perturbation due to the Sun's quadrupole mass moment J2 = 2.0 X 10 — 7 over 
At = 5 yr. The same initial conditions (J2000.0) have been used for both the integrations. The 
state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. The 
integrations have been performed in the ICRF/J2000.0 reference frame, with the mean equinox of 
the reference epoch and the reference {xy} plane rotated from the mean ecliptic of the epoch to 
the Sun's equator, centered at the Solar System Barycenter (SSB). 
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Fig. 25. Difference A|p| = \pp\ — ppj in the numerically integrated EMB-Mars ranges with and 
without the perturbation due to the Sun's Lense-Thirring field over At = 5 yr. The same initial 
conditions (J2000.0) have been used for both the integrations. The state vectors at the reference 
epoch have been retrieved from the NASA JPL Horizons system. The integrations have been 
performed in the ICRF/J2000.0 reference frame, with the mean equinox of the reference epoch 
and the reference {xy} plane rotated from the mean ecliptic of the epoch to the Sun's equator, 
centered at the Solar System Barycenter (SSB). 
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one at a 1.6 x 10~ 5 level; moreover, note the similar time evolution of the two 
signals. Concerning J 2 , the Lense-Thirring effect would bias its signal at a 6% level. 
Conversely, if one looks at J2 as a potential source of systematic bias for the recovery 
of the gravitomagnctic effect, the mismodeled signature of the Sun's quadrupolar 
mass moment would be 1.7 times larger than it. An improvement in its knowledge 
by one order of magnitude, as expected from, e.g., BepiColombo, would push its bias 
on the Lense-Thirring signal at 17%. Anyway, it must be noted that their temporal 
evolutions are different. 

5.4. The ring of the minor asteroids and Ceres, Pallas and Vesta 

The range perturbation of the ring of the minor asteroids is reproduced in Figure 
28. Its nominal peak-to-peak amplitude is 40 m; by considering a 30% uncertainty in 
the mass of such a ring, the related mismodeled signal would be as large as 12 m. It 
would impact the recovery of the Schwarzschild and J2 signals at 5 x 10 -5 , 1.7 x 10 _1 
level, respectively, while the Lense-Thirring signature would be swamped. Note that 
the time signature of the minor asteroids is different from the relativistic ones and 
more similar to that due to Ji- 

The nominal perturbation on the range of Mars by Ceres, Pallas, Vesta is shown 
in Figure 27. For previous analytical and numerical investigations of the their impact 
on the motions of the Earth and Mars, see Ref. 42. The peak-to-peak amplitude is 
as large as 1400 m, measurable at a 7 x 10 -5 level by assuming an accuracy of 10 
cm in the ranging device. The aliasing effect of the present-day mismodeling in J2 
and in the unmodeled Lense-Thirring effect is of the order of 5 x 10 -3 , 2.8 x 10~ 3 , 
respectively. On the other hand, it can be shown that the peak-to-peak amplitude 
of the mismodeled signature of the main asteroids considered here is up to 14 m. It 
would overwhelm the Lense-Thirring effect, while the systematic relative uncertainty 
on the Schwarzschild and J2 range perturbations would be 5.6 x 10 _5 ,2 x 10 _1 , 
respectively. 
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Fig. 26. Difference A|p| = \pp\ — |pr| in the numerically integrated EMB-Mars ranges with and 
without the nominal perturbation due to the ring of minor asteroids with 29 m r i ng = 1 X 10 — 10 Mq 
and -R r i n g = 3.14 a.u. over At = 5 yr. The same initial conditions (J2000.0) have been used for 
both the integrations. The state vectors at the reference epoch have been retrieved from the NASA 
JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 reference frame, 
with the ecliptic and mean equinox of the reference epoch, centered at the Solar System Baryccnter 
(SSB). 
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At the end, it must be remarked that the present analysis of the impact of as- 
teroids on Mars should be regarded just as necessarily incomplete and preliminary. 
Indeed, over timescales of the order of its orbital period 1 , the orbit of Mars is no- 
tably affected by a larger number of different minor bodies 43,44,45 . On the contrary, 
appreciable effects of them on Mercury and Venus occur on timescales some decades 
long 46,47 . Modern ephemerides like 30,48 , e.g, INPOP06-INPOP08 include the grav- 
itational influences of up to 300 most perturbing asteroids of the Martian orbit. 
Limiting to a ring model may lead to inaccuracies. Anyway, recent developments 49 
may have somewhat mitigated such a risk, especially over timescales of a few years. 
Accurate determinations of the masses of about hundred minor asteroids 50,51,52,53 
is one of the main goals of GAIA-i 54 . Isolating the gravitational perturbation caused 
by a single asteroid on Mars, being strongly correlated and mixed up with those of 
many other asteroids, is a non trivial task. It is outside the scopes of the present pa- 
per, being a possible subject for further, dedicated analyses by independent teams 
of skilful astronomers. 



5.5. The Trans-Neptunian Objects 

More serious is the effect on the Mars range of the ring of the TNOs: it is shown in 
Figure 28. Its pcak-to-pcak amplitude is about 5 m, compatible with the present- 
day range residuals available. Thus, the TNOs, whose mass should be conservatively 
considered as uncertain at a 100% level, would impact the Schwarzschild signal at 
2 x 10~ 5 level: the signatures are different. The J2 effect would be biased at a 7% 
level, while the Lense-Thirring one would be overwhelmed by the TNOs, although 
their patterns are not equal. 



'It amounts to about 1.9 yr, i.e. it is comparable to the expected time span for PLR operations. 
J It is an astrometric spacecraft-based ESA mission whose launch is scheduled for 2012. See on the 
WEB http: / / sci.esa.int/science-e/www/area/index.cfm?farcaid=26. 




Fig. 27. Difference A|p| = |pp| — \pu\ in the numerically integrated EMB-Mars ranges with 
and without the nominal perturbation due to 32 Ceres, Pallas, Vesta over At = 5 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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5.6. Violation of the Strong Equivalence Principle 

Moving to exotic effects, a violation if SEP driven by rj = 10~ 5 is shown in Figure 29. 
Its peak-to-peak amplitude is about 5 cm. Also in this case, it would be undetectable, 
and it would be totally swamped by the other dynamical effects considered, modeled 
or not. 



5.7. Secular variation of the Newtonian constant of gravitation 

The effect on the Mars range of a variation of G as large as that determined by 
Pitjeva in Rcf. 22 is reproduced in Figure 30. It is as large as 1 m, i.e. it would be 
barely detectable because it would likely be removed from the signal when fitting 
the initial conditions. Moreover, it would be overwhelmed by the aliasing effects of 
all the other competing dynamical forces considered so far. Suffices it to say that 
the TNOs signal is 5 times larger; the present-day uncertainties in the masses of 
Ceres, Pallas and Vesta yield a signal of up to 14 m. 




Fig. 28. Difference A|p| = \pp\ — ppj in the numerically integrated EMB-Mars ranges with and 
without the nominal perturbation due to the ring of Trans-Neptunian Objects with 22 m r ; ng = 
5.26 x 1O _8 M and i? ring = 43 a.u. over At = 5 yr. The same initial conditions (J2000.0) have 
been used for both the integrations. The state vectors at the reference epoch have been retrieved 
from the NASA JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 
reference frame, with the ecliptic and mean equinox of the reference epoch, centered at the Solar 
System Barycenter (SSB). 




Fig. 29. Difference A|p| = |pp| — ppj in the numerically integrated EMB-Mars ranges with and 
without the nominal perturbation due a violation of SEP for rj = 10 — 5 over At = 5 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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5.8. The Pioneer Anomaly 

Figure 31 depicts the indirect effect of the Pioneer anomaly on the Mars range. Its 
peak-to-peak amplitude is 30 cm. Also in this case, such a potential exotic effect 
would likely be too small to be realistically detected even with a future, advanced, 
cm-level ranging system. The systematic bias due to the other standard Newtonian 
and relativistic effects would be largely overwhelming. 

5.9. Planet X 

Finally, we consider the action of planet X in Figure 32-Figure 33. Its peak-to-peak 
amplitude is 10 — 20 m; it is not in contrast with the present-day range residuals 
from the Martian spacecraft, also because X has not been explicitly modeled in 
producing them, and part of its signature could have been removed from the signal 
in fitting the initial conditions. Such a range perturbation could be measured with 
a future cm-level ranging system. The bias due to the TNOs would be 2 — 4 times 
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Fig. 30. Difference A|p| = \pp\ — |ppj in the numerically integrated EMB-Mars ranges with and 
without the nominal perturbation due to a secular variation of G as large as 22 G/G = — 5.9x 10 — 14 
yr -1 over At = 5 yr. The same initial conditions (J2000.0) have been used for both the integrations. 
The state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. 
The integrations have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and 
mean equinox of the reference epoch, centered at the Solar System Barycenter (SSB). 
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Fig. 31. Difference A|p| = \pp\ — |pr| in the numerically integrated EMB-Mars ranges with and 
without the nominal perturbation due to a Pioneer-like constant and unform radial acceleration 
of 37 Api = 8.74 X 10~ 10 m s~ 2 acting upon Uranus, Neptune, Pluto, Eris over At = 5 yr. The 
same initial conditions (J2000.0) have been used for both the integrations. The state vectors at 
the reference epoch have been retrieved from the NASA JPL Horizons system. The integrations 
have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of 
the reference epoch, centered at the Solar System Barycenter (SSB). 
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smaller, and their signature would be different from that of X. Concerning the 
Lense-Thirring effect, if not modeled it would not mimic the action of X, and its 
magnitude would be 2.5 — 5 times smaller that that of X. The mismodeled signal 
due to the Sun's J2 would be about 1.4 — 3 times smaller, but the S/N ratio would 
become more favorable after the expected improvements by one order of magnitude 
in our knowledge of Ji . 



6. Earth-Jupiter range 

Concerning the Earth- Jupiter range, Figure B-14, pag. 18 of Ref. 21 tells us that 
there are just a few sparse points from the encounters with the Pioneer 10/11, Voy- 
ager 1/2, Ulysses and Cassini probes whose accuracy is of a few km, apart from the 
Ulysses' point. 
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Fig. 32. Difference A|p| = pp — ppj in the numerically integrated EMB-Mars ranges with and 
without the perturbation due to hypothetical remote planet X lying almost in the ecliptic with 
minimum tidal parameter 40 ICx = 1-6 X 10 -26 s — 2 over At = 5 yr. The same initial conditions 
(J2000.0) have been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference epoch, 
centered at the Solar System Barycenter (SSB). 
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Fig. 33. Difference A|p| = \pp\ — |pr| in the numerically integrated EMB-Mars ranges with and 
without the perturbation due to hypothetical remote planet X lying almost in the ecliptic with 
maximum tidal parameter 40 /C x = 2.7 X 1CT 26 s" 2 over At = 5 yr. The same initial conditions 
(J2000.0) have been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference epoch, 
centered at the Solar System Barycenter (SSB). 
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Jupiter will be orbited in the next years by the approved Juno spacecraft 55 : the 
nominal mission duration is 1 yr. Europa Jupiter System Mission, or Laplace, 
(Laplace/EJSM) is an unmanned mission k jointly proposed by NASA and ESA 
in the framework of the Cosmic Vision 2015 — 2025 programme for the in-depth 
exploration of jovian moons with a focus on Europa, Ganymede and the Jupiter's 
magnetosphere; a possible launch date would be in 2020. Continuous ranging to 
both the spacecraft would certainly improve our knowledge of the Earth- Jupiter 
distance, although it is not possible to give a precise figure for that. 
Determining accurately the motion of Jupiter may also have consequences on our 
precise knowledge of certain important features of the Cosmic Microwave Back- 
ground (CMB) having implications on our views on dark matter and dark energy. 
Indeed, the pointing of the detectors of the Wilkinson Microwave Anisotropy Probe 
(WMAP) is determined empirically from observations of Jupiter. Some results ob- 
tained in such a way have been recently questioned; see Refs. 56 ' 57 and references 
therein about the ongoing debate. 

6.1. The Schwarzs child field of the Sun 

Figure 34 depicts the perturbation caused by the solar Schwarzschild field on the 
Earth- Jupiter range. The peak-to-peak nominal amplitude is 5 x 10 5 m. 

6.2. The oblateness of the Sun 

The effect of the solar oblateness on the jovian range is shown in Figure 35. The 
peak-to-peak nominal amplitude is of the order of more than 100 m. Given the 
current level of uncertainty in the solar J2 , its mismodeled effect amounts to 2 x 10 -4 
of the general relativistic Schwarzschild signal. Note that the temporal pattern is 
quite different with respect to the Schwarzschild one. 



k See on the WEB http://sci.esa. int/science-e/www/object/index.cfm?fobjectid=44037. 



Earth— Jupiter range: SSB numerical calculation 




1 2 3 4 5 

Kyi 



Fig. 34. Difference A|p| = |pp| — \pu\ m the numerically integrated EMB-Jupiter ranges with 
and without the perturbation due to the Sun's Schwarzschild field over At = 5 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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6.3. The Lense-Thirring effect of the Sun 

Figure 36 shows the effect of the solar gravitomagnetic field on the Earth- Jupiter 
range. Its peak-to-peak nominal amplitude is as large as about 7 m. Its time signa- 
ture is different from that of J®, whose mismodelled effect would be of the same 
order of magnitude, but is similar to that due to the Schwarzschild field. Conversely, 
the Lense-Thirring effect, if not accounted for, would bias the J® and Schwarzschild 
signals at 7% and 10~ 5 level, respectively. 

6.4. The ring of the minor asteroids and Ceres, Pallas and Vesta 

The nominal range perturbation caused by the ring of minor asteroids is in Figure 
37. The peak-to-peak amplitude is 250 m; the mismodeled component would be 
as large as 75 m. It would overwhelm the Lense-Thirring signature; the alias on 
the Schwarzschild and Ji signals would be 1.5 x 10~ 4 and 75%, respectively. Note 
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Fig. 35. Difference A|p| = |pp| — ppj in the numerically integrated EMB- Jupiter ranges with and 
without the nominal perturbation due to the Sun's quadrupole mass moment J2 = 2.0 X 10~ 7 over 
At = 5 yr. The same initial conditions (J2000.0) have been used for both the integrations. The 
state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. The 
integrations have been performed in the ICRF/J2000.0 reference frame, with the mean equinox of 
the reference epoch and the reference {xy} plane rotated from the mean ecliptic of the epoch to 
the Sun's equator, centered at the Solar System Barycenter (SSB). 
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Fig. 36. Difference A|p| = |pp| — ppj in the numerically integrated EMB-Jupiter ranges with and 
without the perturbation due to the Sun's Lense-Thirring field over At = 5 yr. The same initial 
conditions (J2000.0) have been used for both the integrations. The state vectors at the reference 
epoch have been retrieved from the NASA JPL Horizons system. The integrations have been 
performed in the ICRF/J2000.0 reference frame, with the mean equinox of the reference epoch 
and the reference {xy} plane rotated from the mean ecliptic of the epoch to the Sun's equator, 
centered at the Solar System Barycenter (SSB). 
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that the temporal pattern of the asteroidal ring is quite different from that of the 
previous effects considered. 

The effect of Ceres, Pallas and Vesta is illustrated in Figure 38. The nominal peak- 
to-peak amplitude is of the order of 1000 m; given the 3 — 0.6% level of uncertainty 
in the masses of such major asteroids 32 , their mismodeled signature would amount 
to tens m. 

Such effects would bias the general relativistic Schwarzschild signal at a 6 x 10~ 5 
level, being 4 times larger than the Lense-Thirring one. Concerning the solar J®, 
the mismodelled action of the major asteroids would represent up to 30% of its 
signature. However, also in this case, the temporal pattern is different from the 
other ones. 

6.5. The Trans-Neptunian Objects 

Figure 39 is dedicated to the TNOs. It shows a signature with a nominal amplitude 




Fig. 37. Difference A|p| = |pp| — ppj in the numerically integrated EMB- Jupiter ranges with and 
without the nominal perturbation due to the ring of minor asteroids with 29 m r ; n g = 1 X 10 — 10 Mq 
and -R r i n g = 3.14 a.u. over At = 5 yr. The same initial conditions (J2000.0) have been used for 
both the integrations. The state vectors at the reference epoch have been retrieved from the NASA 
JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 reference frame, 
with the ecliptic and mean equinox of the reference epoch, centered at the Solar System Baryccnter 
(SSB). 




Fig. 38. Difference A|p| = |pp| — |ppj in the numerically integrated EMB-Jupiter ranges with 
and without the nominal perturbation due to 32 Ceres, Pallas, Vesta over At = 5 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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of 80 m exhibiting a steadily increasing pattern. In regard to their corrupting action 
on the signals of interest, by assuming a 100% uncertainty they would bias the 
Schwarzschild and J2 effects at a 1.6 x 10~ 4 and 80% level, respectively. 

6.6. Violation of the Strong Equivalence Principle 

In Figure 40 the effect of a SEP violation as large as r\ = 10 -5 on the Earth- 
Jupiter range is plotted. Its peak-to-peak amplitude is 0.2 m, far too small to 
be detected in any foreseeable future also because it would be swamped by the 
competing Einsteinian and Newtonian signals. 

6.7. Secular variation of the Newtonian constant of gravitation 

The effect of a time variation of the Newtonian constant of gravitation as large 
as 22 G/G = —5.9 x 10~ 14 yr -1 is plotted in Figure 41. Its nominal peak-to-peak 
amplitude is about 2 — 2.5 m; also in this case, it is unlikely that any reasonable 
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Fig. 39. Difference A|p| = |pp| — |ppj in the numerically integrated EMB- Jupiter ranges with 
and without the nominal perturbation due to the ring of Trans-Neptunian Objects with 22 m r i ng = 
5.26 x 1O _8 M and H ring = 43 a.u. over At = 5 yr. The same initial conditions (J2000.0) have 
been used for both the integrations. The state vectors at the reference epoch have been retrieved 
from the NASA JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 
reference frame, with the ecliptic and mean equinox of the reference epoch, centered at the Solar 
System Barycenter (SSB). 
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Fig. 40. Difference A|p| = |pp| — |/7r| in the numerically integrated EMB-Jupitcr ranges with and 
without the nominal perturbation due to a violation of SEP according to rj = 10 — 5 over At = 5 yr. 
The same initial conditions (J2000.0) have been used for both the integrations. The state vectors 
at the reference epoch have been retrieved from the NASA JPL Horizons system. The integrations 
have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of 
the reference epoch, centered at the Solar System Barycenter (SSB). 
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improvements in the ranging to Jupiter may allow for a detection of such a putative 
tiny effect. Moreover, it would be largely overwhelmed by the other relativistic and 
classical dynamical signals. 

6.8. The Pioneer Anomaly 

Figure 42 shows the indirect effect of the standard form of the putative Pioneer 
anomaly on the jovian range through its direct action on Uranus, Neptune and 
Pluto. Its nominal effect is as large as 5 m, likely undetectable. Also in this case, 
the other standard dynamical effects would be larger. 

6.9. Planet X 

The impact of a putative, distant planet X on the range of Jupiter is shown in 
Figure 43-Figure 44. They refer to the minimum and maximum values of the tidal 
parameter K.x of X as derived from the anomalous perihelion precession of Saturn 
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Fig. 41. Difference Ajp| = |pp| — ppj in the numerically integrated EMB- Jupiter ranges with and 
without the nominal perturbation due to a secular variation of G as large as 22 G/G = —5.9 X 10~ 14 
yr — 1 over At = 5 yr. The same initial conditions (J2000.0) have been used for both the integrations. 
The state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. 
The integrations have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and 
mean equinox of the reference epoch, centered at the Solar System Barycenter (SSB). 




Fig. 42. Difference A|p| = |pp| — ppj in the numerically integrated EMB-Jupiter ranges with and 
without the nominal perturbation due to a Pioneer-like constant and unform radial acceleration 
of 37 Api = 8.74 X 10 -10 m s -2 acting upon Uranus, Neptune, Pluto, Eris over At = 5 yr. The 
same initial conditions (J2000.0) have been used for both the integrations. The state vectors at 
the reference epoch have been retrieved from the NASA JPL Horizons system. The integrations 
have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of 
the reference epoch, centered at the Solar System Barycenter (SSB). 
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of —6 ± 2 mas cty -1 preliminarily obtained from an initial data processing of the 
Cassini data themselves analyzed with the EPM ephemerides. The peak-to-peak 
amplitudes of the signals of X amount to 120 — 250 m. Such a putative effect is 
roughly comparable to the (nominal) signatures induced by the Sun's quadrupole, 
the minor asteroids and the TNOs. However, the signature of X is different from 
those of such competing effects, apart from the ring of minor asteroids (Figure 37). 



7. Earth-Saturn range 

After the Cassini spacecraft started its "grand tour" of the Saturnian system, it 
has been possible to drastically increase the accuracy of the orbit determination of 
the ringed planet through direct ranging to Cassini itself. Figure B-20 of Ref. 21 
shows the range residuals of Saturn from 2004 to 2006 constructed with the DE421 
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Fig. 43. Difference A|p| = |pp| — ppj in the numerically integrated EMB- Jupiter ranges with and 
without the perturbation due to hypothetical remote planet X lying almost in the ecliptic with 
minimum tidal parameter 40 Kx = 1-5 X 10~ 26 s — 2 over Af = 5 yr. The same initial conditions 
(J2000.0) have been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference epoch, 
centered at the Solar System Barycenter (SSB). 
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Fig. 44. Difference A|p| = |pp| — |ppj in the numerically integrated EMB-Jupiter ranges with and 
without the perturbation due to hypothetical remote planet X lying almost in the ecliptic with 
maximum tidal parameter 40 K% = 2.7 X 10 -26 s -2 over At = 5 yr. The same initial conditions 
(J2000.0) have been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference epoch, 
centered at the Solar System Barycenter (SSB). 
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cphcmeridcs from Cassini normal points; processing of extended data records of 
Cassini is currently ongoing, so that we will consider an integration time span of 5 
yr. The range residuals of Figure B-20 in Ref. 21 are accurate at 10 m level. Also a 
pair of range residuals from close encounters with Voyager 1 (1980) and Voyager 2 
(1982) are shown: they arc almost one order of magnitude less accurate. 

7.1. The Schwarzs child field of the Sun 

Figure 45 shows the Schwarzschild perturbation of the range of Saturn over At = 5 
yr. Its peak-to-peak amplitude amounts to 5 x 10 5 m. This implies a measurement 
accuracy of 2 x 10 -5 , given the present-day level of uncertainty in the Cassini ranging 
residuals. 

7.2. The oblateness of the Sun 

The Sun's oblateness effect on the Saturn range, computed in a solar equatorial 
frame, is depicted in Figure 46. Its nominal peak-to-peak amplitude is about 100 
m. Thus, the relative accuracy in measuring it is somewhat modest, amounting to 
just 1 x 10 _1 . A 10% uncertainty in the solar J2 implies a mismodeled signal of 10 
m. It represents a bias of 2 x 10 -5 on the Schwarzschild signature, which, however, 
has a different pattern. 

7.3. The Lense-Thirring effect of the Sun 

The Lense-Thirring range perturbation for Saturn, computed in a frame with the 
z axis aligned with the Sun's spin axis, is illustrated in Figure 47. Its peak-to- 
peak amplitude is approximately 7 m, too small to be detected with the present- 
day ranging accuracy to Cassini. It would affect the Schwarzschild perturbation at 
1 x 10 -5 level, while the J2 effect would be biased at 7 x 10~ 2 level. The time 
signature of the gravitomagnetic shift is similar to the larger relativistic effect, but 
it is different from that due to the Sun's oblateness. 
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Fig. 45. Difference A|p| = |pp| — \pn\ in the numerically integrated EMB-Saturn ranges with 
and without the perturbation due to the Sun's Schwarzschild field over At = 5 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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7.4. The ring of the minor asteroids and Ceres, Pallas and Vesta 

The nominal effect of the ring of the minor asteroids on the Saturnian range is in 
Figure 48. Its peak-to-peak amplitude is as large as 80 m, barely detectable with 
the current Cassini ranging. The corresponding mismodeled signal would amount 
to 24 m, compatible with the currently available Cassini range residuals spanning 
2 yr. Such a source of systematic alias impacts the Schwarzschild perturbation 
at a 5 x 10~ 5 level, while the J2 signal is biased by it at a 2 x 10 _1 level. The 
Lense-Thirring signal is overwhelmed by the alias due to the minor asteroids. The 
nominal perturbation on the Saturn range by the combined action of Ceres, Pallas, 
Vesta is computed illustrated in Figure 49. The peak-to-peak amplitude amounts 
to 1750 m, measurable at about 6 x 10~ 3 according to the current level of accuracy 
of the ranging to Cassini. On the other hand, the peak-to-peak amplitude of the 
mismodeled signature of the three asteroids, whose dynamical action is included in 
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Fig. 46. Difference A|p| = |pp| — ppj in the numerically integrated EMB-Saturn ranges with and 
without the nominai perturbation due to the Sun's quadrupole mass moment J2 = 2.0 X 10 — 7 over 
At = 5 yr. The same initiai conditions (J2000.0) have been used for both the integrations. The 
state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. The 
integrations have been performed in the ICRF/J2000.0 reference frame, with the mean equinox of 
the reference epoch and the reference {xy} plane rotated from the mean ecliptic of the epoch to 
the Sun's equator, centered at the Solar System Barycenter (SSB). 
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Fig. 47. Difference A|p| = |pp| — ppj in the numerically integrated EMB-Saturn ranges with and 
without the perturbation due to the Sun's Lense-Thirring field over At = 5 yr. The same initial 
conditions (J2000.0) have been used for both the integrations. The state vectors at the reference 
epoch have been retrieved from the NASA JPL Horizons system. The integrations have been 
performed in the ICRF/J2000.0 reference frame, with the mean equinox of the reference epoch 
and the reference {xy} plane rotated from the mean ecliptic of the epoch to the Sun's equator, 
centered at the Solar System Barycenter (SSB). 
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the dynamical models of all the modern ephemerides, is just 6 m. It would overwhelm 
the Lense-Thirring effect, while the relative bias on the Schwarzschild and Ji range 
perturbations would be 4 x 1CP 5 ,2 x 10 _1 , respectively. 

7.5. The Trans- Neptunian Objects 

The action of the TNOs on the Saturn range is more effective. Figure 50 shows 
that it amounts to 200 m. In principle, it would be measurable by the Cassini 
ranging. The absence of such a signal in the present-day range residuals may be 
due to the fact that the dynamical action of the TNOs was not modeled in the 
DE421 ephemerides used to produce them, so that it is likely that part of the 
TNOs signature, amounting to 50 m over 2 yr, has been removed due to the fitting 
of the initial conditions. The Schwarzschild effect is aliased by them at a 4 x 10~ 4 
level, while the J2 and the Lense-Thirring signals would be swamped by the TNOs. 
It is interesting to note that their temporal evolution is different from that of the 
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Fig. 48. Difference A|p| = |pp| — ppj in the numericaiiy integrated EMB-Saturn ranges with and 
without the nominal perturbation due to the ring of minor asteroids with 29 m,;^ = 1 X 10 — 10 Mq 
and -R r i n g = 3.14 a.u. over At = 5 yr. The same initial conditions (J2000.0) have been used for 
both the integrations. The state vectors at the reference epoch have been retrieved from the NASA 
JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 reference frame, 
with the ecliptic and mean equinox of the reference epoch, centered at the Solar System Baryccntcr 
(SSB). 
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Fig. 49. Difference A|p| = |pp| — \pn\ in the numerically integrated EMB-Saturn ranges with 
and without the nominal perturbation due to 32 Ceres, Pallas, Vesta over At = 5 yr. The same 
initial conditions (J2000.0) have been used for both the integrations. The state vectors at the 
reference epoch have been retrieved from the NASA JPL Horizons system. The integrations have 
been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the 
reference epoch, centered at the Solar System Barycenter (SSB). 
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minor asteroids. 

7.6. Violation of the Strong Equivalence Principle 

Figure 51 shows the effect of a violation of SEP for r) = 10 -5 on the Saturn range. 
Also in this case it is quite negligible because it is as large as 9 cm, almost two 
orders of magnitude smaller that the present-day level of accuracy in the Cassini 
ranging. Moreover, also the other dynamical effects considered would completely 
swamp such a signal. 

7.7. Secular variation of the Newtonian constant of gravitation 

A secular variation of the Newtonian constant of gravitation as large as that reported 
by Pitjeva in Ref. 22 would produce as signal as in Figure 52. Its peak-to-peak 
amplitude is 2.5 m, too small to be detected with the present-day ranging to Cassini. 
Also all the other standard Newtonian and relativistic effects, modeled or not, would 
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Fig. 50. Difference A|p| = |pp| — |ppj in the numerically integrated EMB-Saturn ranges with and 
without the nominal perturbation due to the ring of TVans- Neptunian Objects with 22 m r i ng = 
5.26 x 1O _8 M0 and fi r ing — 43 a.u. over At = 5 yr. The same initial conditions (J2000.0) have 
been used for both the integrations. The state vectors at the reference epoch have been retrieved 
from the NASA JPL Horizons system. The integrations have been performed in the ICRF/J2000.0 
reference frame, with the ecliptic and mean equinox of the reference epoch, centered at the Solar 
System Baryccnter (SSB). 
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Fig. 51. Difference A|p| = |pp| — |ppj in the numerically integrated EMB-Saturn ranges with and 
without the nominal perturbation due to a violation of SEP according to 77 = 10 -5 over At = 5 yr. 
The same initial conditions (J2000.0) have been used for both the integrations. The state vectors 
at the reference epoch have been retrieved from the NASA JPL Horizons system. The integrations 
have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of 
the reference epoch, centered at the Solar System Barycenter (SSB). 
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be orders of magnitude larger; for example, it must be recalled that the mismodclcd 
action of Ceres, Pallas, Vesta is as large as 17 m. 

7.8. The Pioneer Anomaly 

The indirect effect of the Pioneer anomaly on the range of Saturn is depicted in 
Figure 53. It would be undetectable because it would be as large as 4 m over At = 
5 yr. Moreover, such a signal would be totally overwhelmed by the other dynamical 
effects considered so far, both of Newtonian and rclativistic origin. 

7.9. Planet X 

In Figure 54-Figure 55 the impact of a putative, distant planet X on the range of 
Saturn is depicted. Let us recall that they refer to the minimum and maximum 
values of the tidal parameter /Cx of X as derived from the anomalous perihelion 
precession of Saturn of — 6±2 mas cty" 1 obtained from the Cassini data themselves 
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Fig. 52. Difference A|p| = |pp| — |ppj in the numerically integrated EMB-Saturn ranges with and 
without the nominal perturbation due to a secular variation of G as large as 22 G/G = —5.9 x 10 -14 
yr — 1 over At = 5 yr. The same initial conditions (J2000.0) have been used for both the integrations. 
The state vectors at the reference epoch have been retrieved from the NASA JPL Horizons system. 
The integrations have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and 
mean equinox of the reference epoch, centered at the Solar System Barycenter (SSB). 
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Fig. 53. Difference A\p\ = |pp| — ppj in the numerically integrated EMB-Saturn ranges with and 
without the nominal perturbation due to a Pioneer-like constant and unform radial acceleration 
of 37 Api = 8.74 X 10~ 10 m s -2 acting upon Uranus, Neptune, Pluto, Eris over At = 5 yr. The 
same initial conditions (J2000.0) have been used for both the integrations. The state vectors at 
the reference epoch have been retrieved from the NASA JPL Horizons system. The integrations 
have been performed in the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of 
the reference epoch, centered at the Solar System Barycenter (SSB). 





January 19, 2011 1:37 WSPC/INSTRUCTION FILE SSB'multieffects 



38 L. Iorio 



analyzed with the EPM ephemerides. The peak-to-peak amplitudes of the signals 
of X , over At = 5 yr, amount to 300 — 600 m. If we consider just 2 yr, the 
magnitude of the X signatures practically coincide with the present-day accuracy 
of the Cassini ranging. It is interesting to note that they are not in contrast with 
the Cassini range residuals obtained with the DE421 ephemerides of Figure B- 
20 of Rcf. 21. The dynamical action of X was not included in the force models 
of the DE421 ephemerides; as explained before, fitting the initial conditions may 
partially or totally remove an unmodeled signature. The action of the TNOs would 
be a competitor of X, at least on such a timescales. Anyway, the patterns of the 
two signals are different. Instead, J2, given the present-day level of uncertainty 
in it, would not alias the X signature. The Lense-Thirring effect is two orders of 
magnitude smaller. 
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Fig. 54. Difference A|p| = |pp| — ppj in the numerically integrated EMB-Saturn ranges with and 
without the perturbation due to hypothetical remote planet X lying almost in the ecliptic with 
minimum tidal parameter 40 K-x = 1.5 X 10~ 2B s~ 2 over At = 5 yr. The same initial conditions 
(J2000.0) have been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference epoch, 
centered at the Solar System Barycenter (SSB). 




Fig. 55. Difference A\p\ = |pp| — ppj in the numerically integrated EMB-Saturn ranges with and 
without the perturbation due to hypothetical remote planet X lying almost in the ecliptic with 
maximum tidal parameter 40 K.% = 2.7 X 10 -26 s~ 2 over At = 5 yr. The same initial conditions 
(J2000.0) have been used for both the integrations. The state vectors at the reference epoch have 
been retrieved from the NASA JPL Horizons system. The integrations have been performed in 
the ICRF/J2000.0 reference frame, with the ecliptic and mean equinox of the reference epoch, 
centered at the Solar System Barycenter (SSB). 
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8. Summary and conclusions 

In view of a possible future implementation of some interplanetary laser ranging 
facilities accurate to cm-level, we have numerically investigated how the ranges be- 
tween the Earth and all the inner planets plus Jupiter and Saturn are affected 
by certain Newtonian and non-Newtonian dynamical effects by simultaneously in- 
tegrating the equations of motion of all the major bodies of the solar system plus 
some minor bodies of it (Ceres, Pallas, Vesta, Pluto, Eris) in the SSB reference 
frame over a time span 2 years long, apart from Mars, Jupiter and Saturn for which 
we adopted 5 years. One of the major goals of the forthcoming or planned in- 
terplanetary ranging missions like, e.g., BcpiColombo is the accurate (of the order 
of, or better than 10 -6 ) determination of the PPN parameters 7 and (3 discrimi- 
nating various metric theories of gravity. To this aim, it must be recalled that the 
observable used in actual interplanetary ranging tests of post-Newtonian gravity 
consists of two gravitationally affected parts: the one, purely relativistic, connected 
with the Shapiro delay of the propagation of the electromagnetic waves induced by 
the Schwarzschild field of the Sun, and the other one due to the reciprocal Earth- 
planet/spacecraft orbital dynamics; reaching exquisite accuracies in only measuring 
the Shapiro delay is useless if the orbital component is known less accurately. That is 
why we paid attention to several Newtonian perturbations on the planetary ranges 
which may be viewed as sources of systematic uncertainty in the main general rel- 
ativistic Schwarzschild signals of interest. The same holds for other Newtonian and 
non- Newtonian target effects as well. 

It turns out that the general relativistic gravitomagnetic Lense-Thirring effect 
of the Sun, not modeled so far either in the planetary ephemerides or in the analy- 
ses of some spacecraft-based future missions like, e.g., BcpiColombo, does actually 
fall within the measurability domain of future cm-level ranging devices. The more 
favorable situation occurs for Mercury because the relative measurement accuracy 
is of the order of 2 — 5 x 10 -3 by assuming a 4.5 — 10 cm uncertainty in the Earth- 
Mercury ranging, as expected for BcpiColombo over some years of operations. It 
is 2 — 5 x 10 -2 for Venus and 1.2 — 2.5 x 10~ 2 for Mars by assuming the same 
level of uncertainty in the corresponding planetary ranging over 5 yr. In the case 
of Saturn, the peak-to-pcak amplitude of the gravitomagnetic range signal is 7 m 
over 5 yr, too small by about a factor 1.5 with respect to the present-day level 
of accuracy in the ranging to Cassini. If not properly modeled and solved-for, the 
Lense-Thirring effect may also impact the determination of other Newtonian and 
post-Newtonian parameters at a nonnegligible level, given the high accuracy with 
which their measurement is pursued. For example, in the case of BcpiColombo the 
expected accuracy in determining 7 and f3 from the range perturbation due to the 
Schwarzschild field of the Sun is of the order of 10~ 6 ; the Lense-Thirring range per- 
turbation would impact the Schwarzschild one at 4 x 10~ 5 level. Another goal of the 
BcpiColombo mission is a measurement of the Sun's quadrupole mass moment J2 
accurate to 10~ 2 ; the unmodelcd Lense-Thirring effect would bias it at 10 _1 level. 
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From the point of view of a measurement of the Sun's gravitomagnctic field itself, 
it results that a major concern would be the solar oblateness; it should be known 
at a 10 -2 level of accuracy-which is just the goal of BcpiColombo-to allow for a 
reduction of its aliasing impact on the Lensc-Thirring signal down to just 17%. The 
ring of the minor asteroids should be taken into account as well because its mis- 
modeling would impact the gravitomagnctic signal at about 7 x 10 -2 . The lingering 
uncertainty in the masses of Ceres, Pallas, Vesta translates into a potential bias of 
about 3 x 10 -2 . The TNOs, not modeled so far apart from the EPM ephemerides, 
would nominally affect it at a 4.5 x 10 -2 level; it must be considered that there 
is currently a high uncertainty in their mass. However, it must be noted that the 
patterns of such sources of systematic bias are different with respect to the grav- 
itomagnctic one. About Venus and Mars, the measurement of their Lensc-Thirring 
range perturbations would be made difficult by the mismodelcd signals due to J2, 
the ring of the minor asteroids, Ceres, Pallas, Vesta and the TNOs because their 
magnitudes are often as large as, or even larger than, the gravitomagnctic ones, 
although their temporal signatures would be different. 

A Newtonian dynamical effect that has been investigated as a potential source 
of systematic uncertainty in the planetary range signals of interest is the ring of 
the minor asteroids. Its nominal signatures would be detectable because they are 
of the order of 4 m (Mercury), 3 m (Venus), 40 m (Mars), 250 m (Jupiter) 80 m 
(Saturn). However, it is currently modeled in the present-day ephemerides, and the 
uncertainty in its mass is of the order of 30%. Thus, the peak-to-peak amplitudes 
of the mismodelcd effects arc 1.2 m (Mercury), 90 cm (Venus), 12 m (Mars), 75 m 
(Jupiter) and 24 m (Saturn). Concerning BepiColombo and Mercury, the impact of 
such aliasing signals on the Schwarzschild, J2 and Lense-Thirring range perturba- 
tions is 3 x 10 -6 , 4 x 10~ 3 , 7 x 10 -2 , respectively; it is not negligible with respect to 
the expected levels of accuracy. Moreover, it must be recalled that in some proposed 
spacecraft-based tests of post-Newtonian gravity the target accuracies in measuring 
7 and (3 may be as high as 10~ 8 — 10~ 9 . In the case of Venus, the uncertainties in the 
masses of Ceres, Pallas, Vesta translate into a relative bias on the Schwarzschild, 
J2 and Lense-Thirring range perturbations of 4 x 10~ 7 , 1 x 10~ 3 , 3 x 10~ 2 , respec- 
tively. The situation for Mars, which is another possible candidate for implement- 
ing accurate interplanetary ranging, is less favorable because the relative uncer- 
tainties in the three signals of interest are 4.3 x 10~ 5 ,2 x 10 _1 ,3, respectively. 
Saturn should be considered as well in view of possible Cassini ranging tests of 
post-Newtonian gravity. The impact of the mismodeling in the ring of the minor 
asteroids is 5 x 10 -5 for the Schwarzschild range perturbation. The dynamical ac- 
tion of the three major asteroids, i.e. Ceres, Pallas, Vesta, is currently modeled 
in all the modern planetary ephemerides, but the present-day 10~ 3 — 10~ 2 level 
of uncertainty in their masses would induce mismodelcd signatures which cannot 
be neglected with respect to the goal of accurate measurements of the Newtonian 
post-Newtonian parameters of interest. Indeed, in the case of Mercury the system- 
atic uncertainties in the Schwarzschild, J2 and Lense-Thirring range perturbations 
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are 1 x 10~ 6 , 2 x 10~ 3 , 3 x 10 -2 , respectively. For Venus their impact on the first two 
effects is 8 x 10~ 6 , 2 x 10 -2 , respectively, while the gravitomagnctic signature would 
be overwhelmed. It is worse with Mars, for which the alias in the Schwarzschild 
and J2 range perturbations due to the mismodcling in the three large asteroids is 
5 x 10~ 5 ,2 x 10 _1 , respectively. However, it must be recalled that more asteroids 
do actually perturb the orbit of Mars at a nonnegligible level. The gravitoelectric 
signal in the range of Saturn would be affected at a 3.5 x 10 -5 level by them. 

A classical dynamical effect not considered so far in some ephemerides and mis- 
sion analyses is the action of the TNOs, which we modeled as a massive ring. It turns 
out that it may impact some high-precision tests of Newtonian and post-Newtonian 
gravity just at the level of desired accuracy. In the case of BepiColombo, the bias 
of the TNOs on the Schwarzschild and J2 range signals amounts to 2 x 1CP 6 and 
3 x 10~ 3 , respectively. The TNOs' impact on the gravitoelectric and J2 range of 
Venus is as large as 4 x 10 -6 , 1.2 x 10~ 2 , respectively. More effective is their action 
on the Mars range. Indeed, in this case their bias on the Schwarzschild and J2 ranges 
is 2 x 10~ 5 and 8 x 10~ 2 , respectively. It maybe interesting to note that the TNOs 
impact the Schwarzschild range signal of Saturn at 4 x 1CP 4 level; it should be taken 
into account in possible, future Cassini ranging-based tests of post-Newtonian grav- 
ity at the ringed planet. A similar level of bias by the TNOs on the gravitoelectric 
range signal occurs also for Jupiter. 

We also examined other non-Newtonian dynamical effects like a SEP violation 
through the 77 parameter, a secular variation of the Newtonian constant of gravi- 
tation and the indirect effect of the Pioneer anomaly on the inner planets through 
the altered action of the giant planets, putatively acted upon by it, on them. Con- 
cerning the SEP violation, since rj is currently known at 10~ 4 level from LLR, we 
looked at the case r\ = 10~ 5 . It turns out that the largest effect occurs for Mars 
and Saturn amounting to 5 — 9 cm (At = 5 yr); for Mercury and Venus the cor- 
responding SEP signals are as large as 6 mm and 8 mm, respectively, over 2 yr. 
They are realistically too small to be detectable; moreover, they would be com- 
pletely overwhelmed by all the other unmodeled/mismodelcd standard Newtonian 
and relativistic effects. Recently, a statistically significative secular decrease of G 
of the order of 10 -14 yr -1 has been preliminarily reported by E.V. Pitjeva; thus, 
we looked at interplanetary range signals due to such an effect as well. It turns out 
that BepiColombo would, perhaps, be able to detect the corresponding range per- 
turbation for Mercury of 60 cm (peak-to-pcak maximum amplitude) with a relative 
accuracy of 1 — 2 x 10 _1 ; it would impact the Schwarzschild, J2 and Lense-Thirring 
signals at 1 x 10 _6 ,2 x 10 _3 ,3 x 10 -2 , respectively. Conversely, the Lense-Thirring 
and TNOs signals, if not modeled, would severely affect the putative G/G signature; 
the same also holds for the mismodclcd signature of Ceres, Pallas, Vesta. For Venus 
(At = 2 yr) and Mars (At = 5 yr) the range signals would be as large as 0.07 — 1 
m, likely too small to be detectable, also because of the huge biasing action of the 
other competing standard Newtonian and relativistic effects. The G/G range signal 
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for Saturn is 2 m over 5 yr, which is beyond the current capabilities of the ranging 
to Cassini. The indirect effects of the Pioneer anomaly on the interplanetary ranges 
examined are far too small amounting to a few mm for Mercury and Venus (At = 2 
yr) (for Jupiter and Saturn it is up to 4-5 m over 5 yr). 

Finally, we looked at the perturbations on the interplanetary range signals that a 
distant, planetary-sized body X may induce. We considered for its tidal parameter 
the minimum and maximum values obtained by analyzing a putative anomalous 
precession of the perihelion of Saturn recently determined by Pitjcva and Ficnga 
ct al. from preliminary analyses of some years of Cassini normal points with the 
latest versions of the EPM and INPOP ephemerides. Future analyses of extended 
data sets of Cassini should shed more light on the genuine existence of such a 
phenomenon. Our investigations show that the range signals caused by X would 
be well measurable with future cm-level ranging devices. Indeed, the peak-to-peak 
amplitudes of its signatures for Mercury and Venus are 1.5 — 3 m and 3 — 5 m, 
respectively. For Mars, Jupiter and Saturn the peak-to-peak amplitudes of the X 
range perturbations are 10 — 25 m, 120 — 250 m and 300 — 600 m, respectively, over 
At = 5 yr. The figures for Saturn are quite larger than the present-day level of 
accuracy in the ranging to Cassini. 

Table 1 summarizes our results. 

Our analysis can also be extended to other putative, nonconventional gravita- 
tional accelerations induced by modified models of gravity. Moreover, it may be 
helpful in analyzing also other scenarios involving ranging to spacecraft not neces- 
sarily orbiting a given planet or satellite. Indeed, the ASTROD project involves the 
use of three spacecraft ranging coherently with one another using lasers: one should 
be located near the Earth at one of the Lagrange points L1/L2, while the other two 
should move along separate solar orbits. The ASTROD I concept relies upon laser 
ranging from laser stations on the Earth to one spacecraft in solar orbit. LATOR is 
based on the use of a laser transceiver terminal on the International Space Station 
(ISS) and two spacecraft placed in ~ 1 a.u. heliocentric orbits. For such missions 
the expected level of accuracy in determining (3 and 7 is of the order of 10~ 8 — 10~ 9 
after some years of operations. 

At the end, let us elucidate certain limitations of the present study and trace 
some possible routes for further investigations. Our analysis pretends to be neither 
complete nor definitive because, for example, it only accounts for the actions of 
the three major asteroids. This is particularly true for Mars whose orbit is actually 
perturbed by a much larger number of such small bodies; it is well known that mod- 
ern ephemerides, built up with larger human, material and computational resources 
than those available for the present study, include up to 300 biggest asteroids acting 
in a nonnegligible way on the red planet. Another issue necessarily left out by our 
analysis is a complete investigation of the effective detectability of the dynamical 
effects considered. To this aim, it would have been necessary to implement numeri- 
cal integrations fitting the initial conditions to the real observations as well since in 
the actual data processing such a procedure may absorb, to an extent to be quanti- 



January 19, 2011 1:37 WSPC/INSTRUCTION FILE SSB'multieffects 



Effects of Standard and Modified Gravity on Interplanetary Ranges 43 

tatively determined, the signature one is interested in. In this framework, also the 
impact of the noise in the observations on the measurability of the different dy- 
namical effects should be addressed, along with their mutual separability through 
a covariance analysis. However, such important tasks, which are outside the scopes 
of the present analysis, may be the goal for further work by skilful independent 
teams of astronomers routinely engaged in producing even more and more accurate 
ephemerides. 

Table 1. Maximum pcak-to-pcak nominal amplitudes, in m, of the Earth-planet range 
signals due to the dynamical effects listed in the left column for Mercury, Venus, 
Mars, Jupiter and Saturn. The time span used for Mercury and Venus is At = 2 yr, 
while for Mars, Jupiter and Saturn it is At = 5 yr. We adopted the standard value 29 
J2 = 2.0 X 10~ 7 for the quadrupole mass moment of the Sun, while for its angular momen- 
tum we used 27,28 S = 190.0 X 10 39 kg m 2 s~ 1 from helioseismology. For the ring of the 
minor asteroids we used 29 m ring = 1 x lO~ lo M ,R ri ng = 3.14 a.u., while for the TNOs, 
modeled as massive ring as well, we adopted 22 m r i ng = 5.26 X 10 -8 Mq, i? r ing = 43 
a.u. The masses of the major asteroids Ceres Pallas, Vesta have been retrieved from 
Ref. 32. The magnitude of the SEP violation investigated is n = 10 -5 . The secular vari- 
ation of G as been accounted for according to Ref. 22 G/G = —5.9 X 10~ 14 yr -1 . The 
Pioneer anomaly has been included in the forces acting on the outer planets with its 
standard magnitude of 37 8.74 X 10 -10 m s — 2 , while for the tidal parameter of X we 
used 40 GM x /r| = (2.1 ± 0.6) X 10" 26 s~ 2 , obtained from the perihelion precession of 
Saturn 29 - 41 . 
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